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Properties like the excitation energy with respect to the ground state, moments of inertia, B(E2) 
transition probabilities and stability against quadrupole fluctuations at low spin of the predicted 
superdeformed band of 32 S are studied with the Gogny force D1S using the angular momentum 
projected generator coordinate method for the axially symmetric quadrupole moment. The Self 
Consistent Cranking method is also used to describe the superdeformed rotational band. In addition, 
properties of some collective normal deformed states are discussed. 
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I. INTRODUCTION. 

Recently, a number of papers have addressed the theoretical study of the predicted super deformed (SD) configu- 
ration in the nucleus 32 S by using the mean field approximation at high spin with several flavors of the Skyrme |QJ|] 
and also the Gogny interaction fl3|]. The interest to study the SD configuration in 32 S comes from the fact that this 
SD configuration is thought to be an intermediate case between the strongly deformed cluster structures in very light 
nuclei and the known SD states in the A — 60 region Hg]. In fact, it is interesting to understand the relationship 
between the predicted SD band in 32 S and the 16 O + O quasimolecular rotational states observed in this nucleus 
HQ. On the other hand, many states up to an excitation energy of around lOMeV are known experimentally in this 
nucleus S]. Those states can be interpreted in terms of the shell model with active particles in the sd shell Q and 
also in terms of the algebraic cluster model (l(J. It turns out that some of these states can be interpreted in terms 
of deformed intrinsic configurations and, therefore, they can be used as a test ground to assure the reliability of any 
interaction meant to describe the SD band in this nucleus at the mean field level. 

The purpose of this paper is to study, using the Gogny interaction Jll]] with the D1S parameterization jy), the 
properties of the superdeformed band of the nucleus 32 S focusing on the stability of the superdeformed minimum 
at low spin against quadrupole fluctuations. The reason is that in the theoretical studies mentioned in the above 
paragraph and also in previous mean field calculations with the Skyrme JD| and Gogny forces the SD minimum 
observed in the energy landscape was very shallow rising serious doubts about its ability to hold states at low angular 
momentum when fluctuations in the quadrupole degree of freedom are taken into consideration (see J|] for a discussion 
of this issue) . Obviously, at higher spins the rotational energy makes the SD minimum deeper and therefore much 
more stable against quadrupole fluctuations. As a side product of our calculations and also to check the suitability 
of the Gogny force for this nucleus we have studied the properties of low lying normal deformed states and compared 
them with the available experimental data. To perform the theoretical analysis, we have used the Angular Momentum 
Projected Generator Coordinate Method (AMP-GCM) with the axial quadrupole moment as generating coordinate 
and restricted ourselves to K — configurations. The method allows to obtain an accurate estimate of the excitation 
energy of the superdeformed + state with respect to the ground state. The properties of the superdeformed band 
obtained with the AMP-GCM are also compared to those of a Self Consistent Cranking calculation. The choice of 
the Gogny interaction for this calculation is backed up not only by the results we obtain for low lying excited states 
but also by previous calculations in the context of the Bohr hamiltonian in the /3 and 7 collective variables at zero 
spin for normal deformed states that were used to describe neutron and proton pair transfer reactions [jl5| and proton 
scattering ||l6| data with great success. 

II. THEORETICAL FRAMEWORK. 

In the framework of the Angular Momentum Projected Generator Coordinate Method (AMP-GCM) we have used 
the following ansatz for the K = wave functions of the system 

K) = f dq 20 fi(q 20 )P^ |^feo)) . (1) 
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In this expression \(f(q2o)) is the set of axially symmetric (i.e. K = 0) Hartree-Fock-Bogoliubov (HFB) wave functions 
generated with the constraint {(fi{q 2 o)\ z 2 — l/2(x 2 + y 2 ) \(fi(q 2 o)) = Q20 on the mass quadrupole moment. The HFB 
wave functions have been expanded in an axially symmetric Harmonic Oscillator (HO) basis with 10 major shells 
(220 HO states). The two body kinetic energy correction has been fully taken into account in the variational process. 
The Coulomb exchange part of the interaction has not been included into the variational process but added, in a 
perturbative fashion, at the end of the calculation. Reflection symmetry has been used as a selfconsistent symmetry 
in our HFB wave functions. This is not a real limitation as octupole instability is expected [[| at much higher spins 
than the ones considered in this work. 
The operator 

Ho = J dn4 (P)e-^e-^e-^ (2) 

appearing in Eq. (|l|) is the usual angular momentum projector with the K — restriction jl^] and fi(q 2 o) are the 
"collective wave functions" solution of the Hill- Wheeler (HW) equation 

J dq' 20 H I (q 2 o,q 2 o)fU<l2o) = E l J dq 20 M I (q 20 , g 2 o)/i(<?2o)- ( 3 ) 

In the HW equation we have introduced the projected norm JV I (q 20 , q^) = (f(q 2 o)\ Pqq vfeo))' an ^ the projected 
hamiltonian kernel H 1 (q 2 o,q 2a ) = (<p(q2o)\ HP$ \<p(q 20 )). 

The solution of the HW equation for each value of the angular momentum / determines not only the ground state 
(a = 1), which is a member of the Yrast band, but also excited states (cr = 2,3, . . .) that, in the present context, 
may correspond to states with different deformation than the ground state and/or quadrupole vibrational excitations. 
In order to solve the HW equation it is usually convenient to work in an orthogonal basis given by the states 
j/c 7 ) = (nl)^ 1 / 2 f dq 20 u I k (q 20 )pQ \ip(q 20 )) denned in terms of the quantities u^feo) and which are eigenvectors 
and eigenvalues, respectively, of the projected norm, i.e. / dq^JV 1 (q 2 o, <7 2 o) u fe (#20) = n k u k(l^o)- The correlated 
wave functions are written in terms of the new basis as — ^2 k Ik 1 ) where the amplitudes have been 

introduced. In terms of these amplitudes it is possible to define "collective" wave functions g a,I (q 2 o) — J2k 9k' 1 u k(l2o) 
whose square, contrary to the f^{q 2 o) amplitudes, has the meaning of a probability. In the solution of the HW equation 
a technical difficulty is encountered: for q 2 g values close to sphericity and 1^0, the projected norms M 1 (q 2 o, q' 20 ) get 
very small and, as a consequence, the evaluation of the hamiltonian kernels for those values of q 2 o, q' 20 and I is prone 
to strong numerical instabilities. The most notorious consequence is that the angular momentum projected (AMP) 
energy E 1 (q 2 o) = ^(qw, q2o)/N I (q 2 o, 920 ) can not be accurately computed for q 2 o close to sphericity and 7^0. For 
this reason, whenever the AMP energies are plotted in the next section the values corresponding to q 2 o near sphericity 
will be omitted. However, this difficulty does not pose any problem for the solution of the HW equation because the 
configurations with very small projected norms only contribute to the orthogonal states with very small values 
of nj, and these states are not taken into account in the solution of the HW equation. Let us also mention that details 
pertaining the evaluation of the hamiltonian kernels for density dependent forces are given in [jl^,^9| . 

Finally, it has to be said that one of the drawbacks of the method is that the intrinsic wave functions are determined 
before the projection, i.e. we are using a Projection After Variation (PAV) method. A better way to proceed would 
be to use Projection Before the Variation (PBV) pcfl but this would lead us to a triaxial projection which, for the 
moment, is extremely costly for the full configuration spaces used with the Gogny force (see |2^] for an implementation 
of PBV in small configuration spaces). According to the results of refs. J2^j2^] the PBV energy for strongly deformed 

systems (computed after several approximations) can be written as EpBv(I) = (H) — % J j^ + - 2 j^f ) where Jy 
and Jtv are the Yoccoz (Y) and Thouless-Valatin (TV) moments of inertia, respectively. For the PAV one has to 
replace the TV moment of inertia by the Y one in the last term of the previous formula. These results mean that 
the rotational energy correction at I — is already well described in the PAV but for 1^0 one has to use the PBV 
method. As discussed in |Q the effect of the PBV can be estimated by carrying out Self Consistent Cranking (SCC) 
calculations. When the results of these calculations are compared to those of the AMP in several Mg isotopes it is 
found that the SCC gamma ray energies are quenched by a factor 0.7 with respect to the AMP ones. As we will see 
in the next section the same quenching factor appears for the SD rotational band when the SCC and the AMP results 
are compared. 
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III. DISCUSSION OF THE RESULTS. 



In Figure 1 we present the results of the HFB calculations used to generate the intrinsic states |y(<?2o))- On 
the left hand side of the figure we show the HFB energy (lower panel) along with the /?4 deformation parameter 
(middle panel) and the particle-particle energy E pp = — l/2Tr(Are*) (upper panel) for protons and neutrons as a 
function of the quadrupole moment (720- The energy curve shows a deformed ground state minimum at (720 = 0.46 
(/?2 = 0.19) which is only ISOKeV deeper than the spherical configuration. A very shallow super deformed (SD) 
minimum at (720 = 1-96 (/?2 = 0.73) is also observed at an excitation energy E^ FB (SD) = 9.8bMeV. To study the 
effect of the finite size of the basis in the HFB energies and in the excitation energy of the SD minimum we have 
carried out calculations including 18 major shells for the HO basisQ (1140 states) for both the normal deformed (ND) 
and SD HFB minima and found that the corresponding energies are shifted downwards by 759KeV and 1071-ftTey 
respectively. As a consequence, the excitation energy of the SD minimum gets reduced by 312KeV (a 4% effect) up to 
the value E^ FB (SD) = 9.54MeV. The hexadecupole deformation parameter (3a is seeing to increase with increasing 
quadrupole moments and reach at the SD minimum the rather high value (3a = 0.33. Concerning the particle-particle 
correlation energies E pp we observe that their values for protons and neutrons are nearly identical and they go to 
zero in both the normal deformed and superdeformed minima. This implies that dynamical pairing effects could be 
relevant for the description of both the ND and SD bands. On the right hand side of the figure we have plotted the 
matter density contour plots (at a density po = 0.08 fm~ 3 ) for several values of 920- Only for 520 values greater or 
equal 3.2 b the matter density distribution resembles the one corresponding to two touching w O spherical nuclei. On 
the other hand, the matter distribution corresponding to the SD minimum ((720 = 26) resembles closely the density 
obtained in the two-center harmonic oscillator model by coalescing (i.e. taking the distance between the centers of the 
two harmonic oscillator potentials equal to zero) the configuration corresponding to two separate 16 O nuclei in their 
ground states (see ^ and references therein). In this case and according the Harvey prescription the resulting 32 S 
nucleus is not formed in its ground state but rather in the excited configuration (0) 4 (1) 12 (2) 12 (3) 4 where four particles 
have been promoted from the N — 2 major shell to the N — 3 one. According to the ideas developed by Rae |24| ] 
relating clustering to the appearance of shell gaps in the single particle spectrum the above configuration correspond 
to a deformed nucleus with an axis ratio of 2 : 1. In fact, the matter distribution of the SD minimum has an axis ratio 
z/x — 1.92, a proton mean square radius of 3.66/m and deformation parameters /?2 = 0.73 and (3a — 0.33. To study 
further the connection between these ideas and our HFB results for the SD minimum we have computed the spherical 
shell occupancies v(nlj) = J2 m (f \ c nijm Cnl i m \ f) f° r the intrinsic SD wave function. The quantities v(nlj) give the 
occupancy (or contents) of the HO orbital nlj in the intrinsic wave function \ip). The positive parity level occupancies 
are 3.168, 1.688 and 0.706 for the £^5/2, d 3 / 2 and 2s 1/ / 2 orbitals respectively whereas for the negative parity levels the 
occupancies are 0.992, 0.282, 0.864 and 0.174 for the I/7/2, 1/5/2, 2p 3 / 2 and 2pi/ 2 orbitals respectively (the quantities 
for proton and neutrons are very similar and therefore only the proton values are given). Therefore we have for the 
SD intrinsic state 11.124 particles in the N — 2 major shell and 4.624 in the N = 3 one in good agreement with the 
Harvey prescription. As a consequence of these occupancies, we get contributions to the quadrupole moment both 
from the positive parity (sd) orbitals and the negative parity (pf) ones. The two contributions turn out to be nearly 
the same for the SD intrinsic wave function. These values of the occupancies also imply that for a proper description 
of the SD configuration in terms of the shell model one needs to consider not only the sd shell but also the complete 
fp shell. Finally, let us mention that the occupancies of the negative parity orbitals just mentioned are very small for 
the ND minimum as expected. 

On the left hand side of Figure 2 we have plotted the AMP energy curves for I = 0, ... , 12?i (full lines) along with 
the HFB energy curve (dashed line) as a function of the quadrupole moment. The AMP I = Oh energy curve shows 
more pronounced ND and SD minima than the HFB one and they are located at quadrupole moments (720 = 0.556 and 
1720 = 2.026 for the ND and SD configurations. The excitation energy of the SD minimum with respect to the ground 
state for I = Oh is E^ MP (SD) = 8.22MeV to be compared with the 9.85MeV obtained in the HFB calculation. Let 
us mention that performing the AMP calculations with 18 shells is extremely time consuming and therefore we will 
just use in this case the 312KeV shift obtained in the HFB to account for the effect of the finite size of the basis 
in the excitation energy of the SD band head. We will also use the 759KeV shift in the energy of the ground state 
to estimate the binding energy in the AMP-GCM calculation. If we take into consideration the 312KeV shift the 
excitation energy of the SD minimum in the AMP case becomes E£ MP (SD) = 7.91MeV to be compared with the 
9.54MeV obtained in the HFB case with 18 shells. We notice that for increasing spins, the superdeformed minimum 
gets more and more pronounced and becomes the ground state at spin I = 8h. This value is four units lower than the 



1 This basis is considered by some authors M as almost indistinguisable from an infinite basis for the nucleus considered here. 
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HFB predictions of |||| . Finally, let us mention that the main effect of considering the PB V moments of inertia (see 
section 2) in the AMP energy curves would be the lowering of the I ^ curves but the 1 = reference curve should 
remain unchanged. Therefore, we do not expect changes in our prediction of the excitation energy of the SD band 
head and the angular momentum for which the SD band becomes the Yrast band coming from the effects of PBV. 
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FIG. 1. On the left hand side the HFB energy (lower panel), the /?4 deformation parameter (middle panel) and the parti- 
cle-particle correlation energies E vp for protons are neutrons (upper panel) are depicted as a function of the quadrupole moment 
520 given in barns. On the right hand side, contour plots of the matter distribution corresponding to a density of 0.08/m -3 
and different quadrupole moments are depicted. 



On the right hand side of Figure 2 we show the energies obtained in the AMP-GCM calculations for the four 
lowest-lying solutions of the HW equation (labeled with the subindex a = 1, ... ,4) and spins from zero up to \2U. 

Each level has been placed at a (720 value corresponding to its average deformation (<7 2 o) <T ' / = f dq20 |<7 <r '' r (<3'2o)| <?20- 
The I = Oh projected energy curve has also been plotted to guide the eye. 
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FIG. 2. In the left hand side panel the HFB energy (dashed line) and the angular momentum projected energies up to 
I = 12h are plotted as a function of the quadrupole deformation q2o measured in barns. See section 2 for an explanation of the 
missing point around q20 = 0b. In the right hand side, the four lowest-lying solutions of the AMP-GCM equation are plotted 
for each spin. In both plots the Coulomb exchange energy has not been added. See text for further details. 



We first observe that the three lowest + states are spherical whereas the fourth one is located inside the SD 
minimum and therefore is the band head of the SD rotational band. As a consequence of the quadrupole mixing 
the excitation energy of the SD + with respect to the ground state increases up to E^ MP ~ GCM (SD) — 8.87 MeV 
(8.56MeV for the 18 shells basis) and we also confirm that the + SD state is stable against quadrupole fluctuations. 
We also clearly observe the SD rotational band (O4 , 2^ , 4g , 6^, 8f , lOf and 12+) which becomes Yrast at I = Sh 
and shows an irregularity at I — 2h that is related to the near degeneracy of this SD state with the 2g ND state. We 
also notice that a bunch of other SD states with angular momentum 10?i and 12?i appear. Let us also note that, when 
the Coulomb exchange energy is taken into account the Of ground state energy becomes — 270.87MeV r (— 271.63MeV 
when the effect of the finite size of the basis is taken into account) in good agreement with the experimental binding 
energy of 271.780MeV ||. 

Concerning the normal deformed states, the 2+, 4f , 6f and 8% states have nearly the same value of and they 
could be the members of a rotational band with moderate deformation. Moreover, the 2% and 4^" states are oblate 
and together with the 0^" state could be the member of another moderate deformation band. Experimentally, there 
are many known levels at low excitation energy Ja] either of positive and negative parity. They can be interpreted 
in terms of the Shell Model within the sd shell j9| or the algebraic cluster model of [jl(|. Some of the levels can be 
bunched together as members of K — bands and could be associated to the two bands obtained in our calculations. 
In table 1 we compare the available experimental data on excitation energies and B(E2) transition probabilities for 
both K = bands with our predictions. Taking into account that in our calculations we only take into account the 
quadrupole degree of freedom and that the Gogny force has not been fitted for this region of the periodic table we 
conclude from table 1 that our results are in reasonable agreement with experiment specially for the lowest lying 
states. This fact gives us some confidence on the predictions we make for the SD band. In this respect, we can also 
mention that we obtain a spectroscopic quadrupole moment for the 2f of — 13.29/m 2 that compares well with the 
experimental value of — 14.9/to 2 . Our results are also in good agreement with previous calculations with the 
Gogny force in the context of the Bohr hamiltonian for the (3 and 7 collective variables Jl^] . 
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TABLE I. Experimental and theoretical values for the excitation energies (in MeV) and B(E2, 1 ^> I — 2) transition proba- 
bilities (in e 2 /m 4 ) of the two K = bands of 32 S. Values marked with an star correspond to theoretical predictions found in 
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Coming back to the SD band we have also carried out SCC-HFB calculations for the SD intrinsic state to estimate 
the effect of PBV in the moment of inertia of the SD band (see sec. 2). 
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FIG. 3. On the left hand side panel the gamma ray energies E-,(I) — E(I) — E(I — 2) in MeV for the SD configuration 
are depicted as a function of spin 1(h) for the AMP, AMP-GCM and SCC calculations. In addition, the results of the 
AMP-GCM quenched by a factor 0.7 (AMP-GCM-Q) are also depicted. In the right hand side panel the theoretical results for 
the B(E2,I — > 7 — 2) transition probabilities in units of e 2 b 2 are depicted as a function of spin for the SCC and AMP-GCM 
calculations. 



In the left panel of Figure 3 we have plotted the gamma ray energies E 1 (I) = E(I) — E(I — 2) as a function of 
/ for the SD band obtained with different approaches. Namely, angular momentum projection of the SD intrinsic 
configuration, the result of the AMP-GCM and the SCC-HFB results. We observe that the effect of the quadrupole 
mixing is very small as the AMP and the AMP-GCM E 1 (I) are nearly identical. However, the SCC-HFB results 
are clearly different from the AMP ones. The reason was mentioned before and has to do with the fact that the 
moments of inertia computed in the framework of projection before (Thouless-Valatin) and after (Yoccoz) variation 
differ considerably. This fact has already been observed in similar calculations carried out in the same framework for 
several Mg isotopes |Ts[ |. In those calculations we noticed that the SCC-HFB moment of inertia turned out to be a 
factor 1.4 bigger than the AMP one in all the nuclei considered. The same happens here as can be seen by comparing 
the SCC-HFB curve with the one labeled AMP-GCM-Q which is obtained by quenching the AMP-GCM gamma ray 
energies by a factor 0.7. In the following we will consider the AMP-GCM-Q results for the SD rotational band as 
our "best" prediction. However, we have to take into account that the pairing correlations in the SD minimum are 
negligible and therefore one can expect that dynamic pairing could play an important role in determining the SCC 
moment of inertia. The static moment of inertia obtained in the AMP-GCM-Q (or the SCC-HFB) calculation is rather 
constant as a function of spin and has an average value of lO^Mey^which correspond to k = h 2 /(2J) = 48.5KeV. 

In the right panel of Figure 3 we show the B(E2,I — > / — 2) transition probabilities along the SD band for the 
AMP-GCM calculation and the SCC-HFB ones computed within the rotational approximation. In both cases, the 
bare proton charge has been used. Both transition probabilities agree well up to / = 6h where the AMP-GCM results 
become bigger. This is related to the fact that in the SCC-HFB the Coriolis anti-stretching effect diminishes the 
intrinsic deformation from fa — 0.74 at I = to 02 = 0.68 at I — 12 whereas in the AMP-GCM the quadrupole 
deformation increases with spin (see the right hand side panel of Fig. 2). The Coriolis anti-stretching effect cannot 
be present in a AMP-GCM calculation because the intrinsic states used there are those computed at zero spin. The 
computed B(E2) are highly collective and correspond to around IQQW.u. for the 2 + — > + transition. 
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IV. CONCLUDING REMARKS. 



As a result of our calculations we can conclude that the predicted superdeformed band in S is stable at low 
spin against quadrupole fluctuations in spite of the shallowness of the HFB SD minimum. The effect of angular 
momentum projection and quadrupole mixing reduces the excitation energy of the SD band head with respect to the 
ground state to 8.8TMeV (8.56MeV if the effect of the finite size basis is considered), to be compared to the HFB 
result of 9.85MeV (9.54MeV for the 18 shells basis). We have also found that the SD band becomes Yrast at angular 
momentum / = 8h in our calculations in contrast with the results of several HFB results where it becomes Yrast at 
/ = I2h. The computed moment of inertia of the SD band is of the order of 10. 3h 2 MeV" 1 and the B(E2) transition 
probabilities along this band are very collective, exceeding lOOt-Kw. The SD band head configuration corresponds 
to the promotion of four particles from the sd to the fp shell. As a consequence, the matter distribution of the SD 
configuration corresponds to two coalescent 16 O nuclei in the sense of the two-center harmonic oscillator shell model 
in good agreement with the Harvey prescription. 
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